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The Schwinger model (QED2) with N flavors of massive fermions on a circle of circumference L, 
or equivalently at finite temperature T, is reduced to a quantum mechanical system of N — 1 degrees 
of freedom. With degenerate fermion masses (m) the chiral condensate develops a cusp singularity 
at 8 — ±7r in the limit L — > 00 or T — > 0, which is removed by a large asymmetry in the fermion 
masses. Physical quantities sensitively depend on the parameter mL or m/T, and the m — > and 
L — > 00 (or T — > 0) limits do not commute. A detailed analysis is given for N = 3. 
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Two-dimensional QED has profound resemblance to four-dimensional QCD, including chiral symmetry breaking, 
confinement, instantons, and 9 vacua JjJ. Defined on a circle [|| |J, the model is mathematically equivalent to its 
finite temperature version 0], and recently the model with general fermion content was studied in connection with 
Z n symmetry breaking ||. A subtle difference between the massless and massive theories has been noted there, for 
■ which the non-commutativity of the massless limit and zero temperature limit is an crucial factor. In lattice gauge 
theory, the commutivity of the two limits m — > and L — > 00 is subtle and important. This is particularly true in 
the issue of the triviality of QED4 at strong coupling |7| . 

In a previous paper || we have shown that massive two-flavor QED2 defined on a circle is reduced to the quantum 
I ■ mechanics of a pendulum. The physics is controlled by the strength of the pendulum potential k, which is given in a 
I/"") I large volume L by k ~ mL(eL) 1 / 2 \ cos^#| where e, m, and 8 are the coupling constant, fermion mass, and vacuum 
ON ■ angle, respectively. It was recognized that L — » 00 and m — > limits do not commute. Here we generalize our analysis 
I to the iV-flavor case. 

The N > 2 Schwinger model is distinctively different from the N = 1 model. With massless fermions the spectrum 

Qh. contains N — 1 massless bosons (mesons), and the chiral condensate vanishes, (ipip) — 0. P] With massive fermions, 
the boson masses and chiral condensate are non-vanishing, but have singular dependence on m and 6 at 9 — ±7r in 
the L — > co limit. Q 

Suppose that the fermion and gauge fields obey anti-periodic and periodic boundary conditions, respectively. The 
model, after a Wick rotation, is equivalent to the model on a line at finite temperature T = L^ 1 . 
5_h ' In the bosonization method 
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where C± is the Klein factor. We refer the reader to Q for details. Anti-Periodicity for the fermions is ensured by a 
physical state condition e 27 " p ± | phys) = | phys). Thus p± takes integer values. 
The Hamiltonian in the Schrodinger picture becomes 
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Qw is the Wilson line phase around the circle: A\ = @w(t)/eL. It couples to p±'s through the chiral anomaly [pi. 
0a = 0+ + 0" and II a is its canonical conjugate. 
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In the absence of fermion masses the Hamiltonian is exactly solvable. The spectrum contains one massive field 

AT- 1/2 Ef=10a With a mass \i — (N/ir^^e, and N — 1 massless fields. The mass term H mass = L dx J2 a m a^ ai>a 
however gives various interactions among the zero and <f> modes. We present an algorithm valid for \m a \ <C e to evaluate 
the effects of H mass . We stress that this by no means implies that H meiss can be treated as a small perturbation. On 
the contrary, it dominates over Hq in the L — ► oo or T — ► limit; its effects are quite non-perturbative. 

As flmass commutes with p\ — pt_ , we restrict ourselves to states with (p° — p1_ ) | phys) = 0. Then a complete set 
of eigenfunctions of H is $C ni '-' n ^ ~ u s [O w + {2nJ2 a n a/N)] e^^il+i-) where u s satisfies \{-d 2 x + x 2 )u s = 
(s + i)w s with Qw = (ne 2 L 2 /Ny/ 4 x. The ground states of H are infinitely degenerate for m = • • • = tin- 

It proves to be much more convenient to work in a coherent state basis $ s (<^i, • • • , (/Jjv-i; given by 

{n,r a } 

Htot induces transitions among ^ s (ip a ', 0). However the effect of transitions in the s index is small § and we restrict 
ourselves to s — states. The vacuum state is written as 

$vac(0) = / d<px--- d<p N -i f(ip a ; 0) <S> ((p a ; 0) . (4) 
Jo 

iJ mass significantly alters the vacuum structure of the (f> a modes as well. Its main effect is that the N — 1 previously 
massless fields now acquire finite masses. The vacuum is defined with respect to these physical fields. 

We write mass eigenstate fields as Xa = U aa 4>a- In this physical space {p°\_ — pt) |phys) = 0, the fermion mass 

operator M a — ijj a h{l + j 5 )ip a is (in the Schrodinger picture) 
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B( z ) = —exp('y + --2 f°° dH , I (5) 



B(fxL)^^ 7r N^[e i ^]. § It is easy to find 



Here N^l ■ ■] indicates normal-ordering with reference to a mass /i. We have made use of N^[e l ^ x ] = 
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\^ s {6';ip' a )\H Taass \^ s {6;^ a )) = -S 2 ^(9' 6) \\ 5 2 „((p' b - <p b ) ^ 2m a A a cos tp a 
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A a = e-^ L J] BM^ 2 , (6) 



where (p. N = 6 - Y<a=i <Pa- 



The equation H tot $ V ac(#) = E$ vac (6) becomes: { - + V N ((p)} f(ip) = where 
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V N (<f) = - Jjy—^ Yl maAa C ° S ■ (7) 

The potential Vn(<p) depends, through A a defined in (^), on [i a and U aa which are to be self-consistently determined 
from the wave function f(ip a ;6). In the ^-vacuum (|j), (M a ) g = —(A a /L) {e~ liPa ) f where the /-average is defined 

by (g{<p)) f = J[dv]g(<p)\m\ 2 '. 

Vn((p) has a similar structure to the potential which appears in the effective chiral Lagrangian of QCD. In Witten's 
formalism (TT) the faS are related to the pseudoscalar mesons themselves. The constraint Yla=i fa — appears when 
the chiral anomaly dominates over the quark masses. In our case the Xa' s represent the boson fields, whereas ip a 's are 
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parameters of the coherent state basis. The similar structure emerges as a result of the pattern of SU(N) x SU(N) 
symmetry breaking. 

To find the boson masses /i a , we denote 

Ra \ 8tt (Re 



-ffmass yields Yl a ^„ [e lUaay ^ Xa ]. Expanding \a in H mass and adding the contribution from the chiral anomaly, one 
finds 

U aa I a ,1 2,2 , 3\ 



HX ass = I dx {^x^-^xi + Oix*)} ■ (9) 



Here C/ QQ diagonalizes the matrix 
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and jU a are the eigenvalues. (||), (0), @, and the diagonalization of ( [Toj ) must be solved simultaneously. 

Suppose that all fermion masses are degenerate: m a = m <C e. In this case R a = R and ( e~' liPa ) * = { e~' lip ) ^. One 
can choose U aa such that Ui a = N^ 1 / 2 , fi 2 = fj 2 + R, and (£ = ■•• = ^i 2 N = R. The potential and boson masses are 
reduced to 
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V N (lfi)= - Kg COS (fig 
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The last relation is analogous to the PCAC relation in QCD [fl2|| . The strength ko and fix the potential Vn(<p), and 
thus control the physical behavior. Since B(z) ~ e 7 ;z/47r for z ^> 1 and 5(0) = 1, Kq — * oo (0) as L — ► oo (0). So 
long as to ^ 0, the location of the minimum of the potential determines the physics at L — ► oo. 
The potential Vjv(</?) has a minimum at 



<p 1 =.-. = <p N = —(tl- 2;r 



2tt 



(12) 



[— 7r < < 7r]; its location jumps discontinuously at — 7r from <y9 a = 7r/iV to ip a = —Tt/N. There is a special feature 
in the N = 2 case: as V2(<p) = — 2kq cos ^0 cos(v?i — ^0), the potential vanishes at = ±n. Its behavior is controlled 
by the single parameter k q \ cos |0| (see g]). 

When kq > 1 (or as £ oo), the wave function /(</?) approaches a delta function around the minimum of the 
potential Vjv(y). Hence lirnz^oo ( cosip a ) f = cos(0/iV). For to//i <C 1 one finds from ( |ll| ) 



1 i a i 2W , N-l 

J- / — \ i / _ (7 \ JV + 1 /TO\ JV + 1 



-(^) e = --^cos-j (-) . (13) 

In the opposite limit k <C 1, the wave function is given by / = (2w)~^ N ~ 1 ^ 2 {l + n J2a=i cos fa H } so that 

/ \ /(1 + cos0)k o foriV = 2 , > 

(COS ^ > ^{; foriV>3. (M) 



For iV > 3 and fiL < 1, 



-(^)e = -^ L ^-e-^ L . (15) 
/j, 7r(iV — 1) /x 
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In the intermediate region, \xL 3> 1 ^> [i^L, 



n(N - 1) /J, \ 4tt 

For N = 2, (|l5|) and ([l6]) must be multiplied by a factor 2 cos 2 |#. 

For general values of /xi and m/fi, we have determined £t2 and (ipip)g in the N = 3 case numerically. Fig. 1 shows 

typical wave functions |/(</?)| 2 at k = 0.1, 10 and 9 = 0, 0.999?r. In fig. 2, (if>tp)/(J, at 6» = and /xi = 10 3 is plotted 

as a function of m//i. In fig. 3, ( ip ip ) / fj, at 9 = is plotted as a function of T//i and m//i. In fig. 4, the 9 dependence 

of (i/ji/j)/fi at m/[i = 10~ 3 is depicted for various [iL. 

Several important observations follow. As is evident from (|l3|), (ipip) in the infinite volume limit has fractional 
power dependence on the fermion mass to. However, if the massless limit m — *■ is taken with a finite L, then k,q 

becomes very small (kq <C 1) and (ipip) is given by either ( |T5| ) or (|l6|), which is linear in to. In this limit the fermion 
mass term in the Hamiltonian can be treated as a small perturbation. On the other hand, the effect of finite fermion 
masses in infinite volume is always non-perturbative. The massless and infinite volume limits do not commute with 
each other. The important parameter is kq- We have juxtaposed a plot for kq in fig. 2 from which one can see that 
the crossover in physical behavior takes place around kq = 0.2. 

Our result can be reinterpreted for the Schwinger model on a line at finite temperature by replacing L by T _1 . 

We see that there is no phase transition at finite temperature. || The condensate (ipip) , which is non- vanishing at 
T = 0, smoothly goes to zero at finite temperature. See fig. 3. 

Thirdly all physical quantities are periodic in 9 with period 2ir. At L = oo, (ipip) g has a cusp at 9 — ±tt, [|l0| while 
at any finite L the cusp disappears as shown in fig. 4. The appearance of the cusp is traced back to the discontinuous 
jump in the location of the minimum of the potential Vn (</?)■ 

So far we have concentrated on cases with a symmetric fermion mass term. For a general fermion masses, the 
evaluation procedure is more involved. One important conclusion can be drawn. In the potential Vn(<p) in (0), the 
coefficients of cos ip a are all different in general. The potential for N = 3 is proportional to 

F(ip) = —q cos ipi — r cos (f2 — cos(# — tpi — ip2) (17) 

where q = ra\A\jm-iA-i and r = raiAijra^A-i. We have investigated the location of the minimum of F(ip) as a 
function of 9 with various parameter values (q, r). In the symmetric case (g, r) = (1, 1), the location of the minimum 
moves from (ipi, (^2) = (— ^tt, — to (3^, ^7r) as 9 varies from — 7r to +7r, and jumps to return to the original point 
(— |7r, — ^7r); see fig. 5. Now we add an asymmetry. Several cases are plotted in fig. 5. One can see that so long as the 
asymmetry is small enough, there is a discontinuous jump at 9 = ±7r. However, a sufficiently large mass asymmetry 
removes this discontinuity. For instance, for (q,r) = (1,0.3), the minimum moves from (^i,^) = (0, — 7r) to (0, +tt), 
hence making a closed loop in the ip±—ip2 plane. We have observed that with sufficiently large asymmetry, the minimum 
at 9 = ±7r is located at either (0, 0), (0, ±7r) or (±7T, 0). These three points are related by S3 transformations. 

We conclude that a sufficiently large asymmetry in the fermion masses removes the cusp singularity at 9 = ±7r in 

(ipip) in the L — > 00 limit. A similar conclusion has been drawn in the QCD context by Creutz. Jl3| 
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Figures: 




FIG. 1. Typical wave functions |/(v)| 2 a t k-o = 0.1 and 10, and 9 = and 0.999-7T 
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FIG. 2. The behavior of = —{ipip)/fx versus m/fi at /iL — 10 3 and 6 — 0. The analytic forms, (13) and (16), are also 
displayed. The crossover occurs at Kq w 0.2. 
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FIG. 3. The chiral condensate (ipip)/fj, as a function of temperature T/fj, (or 1/jJ.L) at 6 — 0. 
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FIG. 4. The 9 dependence of {i>ip). In each case m/fi = 0.001. 
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FIG. 5. The location of the minimum of F(tpi, ifi2) in (17). The points in the graph for each q and r pair run from 9 = (at 
the origin) and end at 9 = n. For 9 < 0, the minimum is located at (<pi, ip 2 )e = —{fi, ^2) -9- For (q, r) — (1.5, 7) the minimum 
returns to the origin at 9 — it. 
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